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ABSTRACT 

We estimate the incidence of multiply-imaged AGNs among the optical counterparts 
of X-ray selected point-like sources in the XXL field. We also derive the expected 
statistical properties of this sample, such as the redshift distribution of the lensed 
sources and of the deflectors that lead to the formation of multiple images, modelling 
the deflectors using both spherical (SIS) and ellipsoidal (SIE) singular isothermal mass 
distributions. We further assume that the XXL survey sample has the same overall 
properties as the smaller XMM-COSMOS sample restricted to the same flux limits 
and taking into account the detection probability of the XXL survey. 

Among the X-ray sources with a flux in the [0.5 — 2] keV band larger than 3.0 x 
10“^® erg cm“^ s“^ and with optical counterparts brighter than an r-band magnitude 
of 25, we expect ~ 20 multiply-imaged sources. Out of these, ^16 should be detected 
if the search is made among the seeing-limited images of the X-ray AGN optical 
counterparts and only one of them should be composed of more than two lensed images. 
Finally, we study the impact of the cosmological model on the expected fraction of 
lensed sources. 
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1 INTRODUCTION 

The XXL surveR carried out by the space-based X-ray ob¬ 
servatory XMM-Newton, spans over ~ 2 x 25 square de¬ 
grees with near 10 ks exposure in each field and is expected 
to lead to the detection of ~ 25000 Active Galactic Nu¬ 
clei (AGNs) down to a limiting flux 1 0~^^ erg cm~^ s~ ^ 
in the [0.5 — 2] keV soft X-ray band (iPierre et al.l l2015ll . 
These X-ray data are complemented by multi-wavelength 
data obtained with the Canada-France-Hawaii Telescope 
Legacy Survey (GFHTLS) and with the Blanco telescope 
(Blanco/South Pole Telescope (SPT) Cosmology Survey, 
BGS) in the (near-)optical u\ g, r, i and z bands, down 
to a limiting AB magnitude of ~ 25. Beside the multi-band 
imaging of the XXL fields, there is a very large on-going ef¬ 
fort to obtain optical spectra of XXL sources, through either 
the matching of existing survey catalogues or dedicated spec¬ 
troscopic surveys. Among these spectroscopic data acquisi¬ 
tion programmes, the VIMOS Public Extragalactic Redshift 
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Survey (VIPERS. lA.Guzzo fc Le Fevrd FoiOf) covers most 
of the northern field, the southern field being covered us¬ 
ing the AAOmega multi-object spectrograph on the Anglo- 
Australian Telescope, an instrum ent used for the Ga laxy and 
mass assembly project (GAMA. [Driver et al.|[200^ . 

The completeness of this multi-wavelength database 
over the entire XXL field provides a unique sample to search 
for multiply-imaged AGNs. We have thus initiated such a 
search among the optical counterparts of point-like sources 
in the soft X-ray band. Beside the scientific interest provided 
by each multiply-imaged source, the goal of this project is 
to construct a statistically clean sample of lensed sources 
that will be used, in combination with samples of multiply- 
imaged sources from other recent surveys, to independently 
constrain the cosmological model. 

The choice of the soft X-ray point-like sources is moti¬ 
vated by the higher sensitivity of XMM-Newton in this band. 
Furthermore, this spectral band should contain a larger frac¬ 
tion of type-I AGNs than the hard X-ray. On average, type-I 
AGNs with a detectable optical counterpart are expected to 
have a higher redshift than type-II AGNs (more absorbed 
in the visible and thus more difficult to detect in the optical 
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at high redshift). As higher redshift sources have a higher 
probability of being lensed, this population is more likely to 
undergo gravitational leasing with the formation of multi¬ 
ple images. The better angular resolution achievable in the 
optical domain will allow to unravel the multiply-imaged 
sources. 

A search for gravitational lenses among the optical 
counterparts of X-ray sources has already been carried out 
for a subset of the XX L field, the XMM- Larqe Scale Struc¬ 
ture field (XMM-LSS. lElviv et al.ll2013li . For this smaller 
held, visual inspection of all optical counterparts has been 
done in order to identify the multiply-imaged AGN can¬ 
didates that are now awaiting spectroscopic conhrmation. 
The search for lensed sources in the larger XXL held is in 
progress. 

In this paper, we present a prospective analysis of the 
lensed AGN population detectable within the XXL survey, 
as well as a study of their expected statistical properties. 
In order to perform this analysis, we have reformulated the 
mathematical formalism to study the statistical aspects of 
gravitational lensing, basing the statistical formalism on the 
observables of the source population. 

This paper is structured as follows. We present the 
mathematical approach in Section [5] Namely, we derive the 
expression allowing to calculate the probability for a source 
to be lensed with the formation of multiple images, mod¬ 
elling the dehector population by means of a spherical mass 
distribution and then taking into account the internal el- 
lipticity of the dehector mass distribution. We explain how 
this expression may be averaged over the entire population 
of sources detected in the survey, thanks to the source joint 
probability density, with which we derive the expression of 
the expected fraction of lensed sources in a survey, as well as 
the expected redshift distribution of the lensed sources and 
that of the dehectors. Our simulations also account for the 
inability of the ground-based GFHT and Blanco telescopes 
to resolve multiple images with too small angular separa¬ 
tions. 

In Section [31 we present the observational constraints 
used to estimate the expected properties of the XXL pop¬ 
ulation in the X-ray and optical domains: these were in¬ 
ferred from the deeper (but smaller) XMM-COSMOS held 
llBrusa et al Finally, in Section [4] we present our re¬ 

sults, i.e. the expected number of multiply-imaged sources 
in the XXL and the XMM-LSS helds, as well as the expected 
statistical properties of these lensed sources. We also inves¬ 
tigate the fraction of lensing events as a function of their 
number of lensed images and we investigate how the frac¬ 
tion of multiply-imaged sources changes as a function of the 
cosmological mass density parameter, Q,m- 

2 MATHEMATICAL FORMALISM 
2.1 Lensing optical depth 

Multiple images due to gravitational lensing occur when 
light rays emitted from a background source are deviated 
towards the observer by a foreground dehector located near 
the line of sight. In our case, the amplihed lensed sources 
have to be above the survey hux limit in both the X-ray 
and the optical domains, and the lensed images have to be 


resolved in the latter. The probability for a source to be 
lensed depends on its redshift, its X-ray hux and r-band 
magnitude. In this section, we derive an expression to calcu¬ 
late the probability for a source to undergo a gravitational 
lensing event as a function of its redshift and its X-ray hux 
only, considering the minimal angular separation resovable 
in the r-band. We will analyse the validity of this simplihca- 
tion in Section [5] by formally including the r-magnitude in 
the calculations. 

Let us consider a source with redshift Zs, with an ob¬ 
served hux fx in the X-ray band, and a lens with a mass 
distribution characterised by a set of parameters M ', located 
at an intermediate position along the line of sight at redshift 
Zd- In the lens plane perpendicular to the line of sight, we 
can dehne an area T, (zs, Zd, fx, M') centred on the source 
projected on the lens plane, called the lensing cross section 
in which the presence of a dehector leads to the detection of 
multiple images by the observer (where the multiple images 
are resolved in the r-band GCD frames). The lensing cross 
section is a function of the redshifts of the source and the 
dehector, the X-ray hux fx of the source and the dehector 
mass distribution parameters M' (some mass distributions 
are more efficient at dehecting light rays and thus have a 
larger lensing cross section). 

The probability dr (za, Zd, fx,M') for this source to be 
multiply imaged due to the presence o f a dehector in the 
redshift range [zd, Zd -f dzd] is given by ([Turner et al.lfl983l 

dr = (1 -I- Zd)^ riD {zd, M') ^ E [zs, Zd, fx,M') dzd, (1) 

where no {zd, M') is the volume density in the comoving ref¬ 
erence frame of dehectors characterised by the mass param¬ 
eters M'. The quantity cdt/dzd represents the inhnitesimal 
light-distance element at redshift Zd per dehector redshift 
unit, which, in a h at expanding FLRW universe, is given by 
(e.g. |Peeblesiri993l ') 

where Q.m is the present-day value of the cosmological mass 
density parameter. 

The envelope of the lensing cross sections at diherent 
dehector redshifts Zd dehnes the lensing volume in which 
the presence of a dehector leads to the detection of multi¬ 
ple lensed images of the background source. The probability 
T {zs, fx, M') for a source to be lensed with the formation 
of multiple images can be calculated by integrating Eq. [T] 
over all values of the dehector redshift Zd, which leads to 

r {zs,fx,M') = 

y (1 + Zdf riD {zd, M') ^ E [zs, Zd, fx, M') dzd- 

0 

The dehnition of t {za, fx, M') in Eq. [3] corresponds to an 
optical depth, which for small values can be assimilated to a 
probability. For this reason we refer to r {zs, fx, M') as the 
source lensing optical depth or lensing probability, without 
distinction. 

In Eq.|3] the comoving density ud izd, M') of dehectors 
assumes one type of dehectors with similar characteristics 
dehned by the parameters M'. Considering the mass distri¬ 
bution to be characterised by means of dehectors with a cen¬ 
tral velocity dispersion a, in the range \a, a -\- da], no {zd, a) 
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may be expressed by means of the Velocity Dispersion Func¬ 
tion (VDF) of galaxies 3>o- (cr, Zd) 


no (zd, cr) = (cr, Zd) da. 


(4) 


To take into account the contribution of all galaxies with dif¬ 
ferent central velocity dispersions, Eq.[3]has to be integrated 
over a. The VDF in Eq. |4]can be either measured directly 
or inferred from the Luminosity Function of the deflecting 
galaxies, using the mean Faber-Jackson or Tully-Fischer re- 
lationshi p, depending on t he type of galaxies considered. 
However. ISheth et al.l ll2003l i have shown that neglecting the 
dispersion of the Faber-Jackson relationship leads to a wrong 
estimate of the VDF. We thus use the VDF de termined di¬ 
rectly through observations (thereby follow ing Shethj^ajJ 
200alMitchell et al.ll2005l. IChoi et aj.ll2007l and lChae|[20l^ . 
Sheth et al. ( 2003h and Choi et al.l 1 200'^ have shown that 
the VDF of early- and late-type galaxies is well fitted by the 
modified Schechter function 


[a, Zd) da = <E>* 


exp 


/3 da 


(5) 

where and a, are the characteristic number density and 
central velocity dispersion, respectively, a and j3 are the 
slope coefficients of the VDF for low and high values of a, 
respectively, and where F (x) is the Gamma function. 

Although less numerous, early-type galaxies are much 
more efficient deflectors than late-type ones, which tend 
to form multiple images with smaller angular separa¬ 
tion. Late-type galaxies contribute by less than 10% to 
the gravitational lensing events in a typical sample of 
lensed AGNs selected i n the optical or near infrared 
i Fukugit a fc Tur neijl99ll. Maoz fc Rixl (ll993l L lKeeton et al.l 
1998li . iKochanek et al.l feoodi L although in radio-selected 


samples, thanks to the better angular resolution of the sur¬ 
vey, the fraction of lensing events formed by late-type deflec¬ 
tors may be higher (e.g. the CLASS survey where at least 
5 of 22 deflectors are late-type galaxies, see iBrowne et all 
l2003l 'l. In the present work, because the multiple images will 
be searched for in the SDSS r-band, we consider the deflector 
population to be only composed of early-type galaxies, and 
we will study in Section 3] the validity of this assumption. 

Constraints from strong-lensing statistics on the evolu¬ 
tion of the VDF of early-type galaxies show very little evo- 
luti on or are consistent with a no-evolution ass umption (see 
e.g.lChaell200l lOfek et al.ll200lIChadl2010l and lOguri et al.l 
l2oiar Consequently, throughout this work we assume the 
deflector VDF to be constant with redshift in the comoving 
reference frame, or in other words, that there is no impact 
of the evolution of the deflector population on the VDF and 
we use the value of the VDF param eters determined in the 
local universe bv IChoi et al.l (l2007l L i.e. 

[-F., a*, a, /3] = [8 X 10"® h® Mpc■^ 161 kms"\ 2.32, 2.67] . 

The lensing cross section in Eq. [T] depends on the de¬ 
flector mass distribution parmeters M'. For a deflector mass 
distribution modelled as a Singular Isothermal Sphere (SIS) 
(a spherical mass distribution with a volume density scaling 
as r~^), the mass distribution is characterised by the line 
of sight velocity dispersion a and th e lensing cross sect ion 
E (za, Zd, fx, a ) can be defined as isee lTurner et al1ll984l for 
definition and I Claeskend 0*9991 for the formalism followed in 


this paper) 


E{zs,Zd,fx,a) 



Nf^ Ux/A{y)) 
Nf^ (fx) 


dy, 


( 6 ) 

where we have introduced the source vector y = (yi,y 2 ) 
which Cartesian coordinates are projected on the deflector 
plane, and normalised to the scale factor bo (i.e. the Einstein 
radius in the deflector plane). A{y) is the total amplifica¬ 
tion of the lensed images formed for a source located at the 
position y, i.e. the sum of the amplification moduli of each 
lensed image and Nf^ (fx) is the differential number counts 
function (DNCF) as a function of the source flux fx in the 
X-ray band . For an SIS de flector, the scale factor is given 
by (see e.g. IClaeskeri^ll999h 


bo {a, Zd, Zs,UM) = d-TT 



DodDds 

Dos 


(7) 


where Dqd, Dds and Dos represent the different angular 
diameter distances between the observer, the deflector and 
the source, c is the speed of light, and UM = (f2m,F7o) is 
a set of parameter values characterising the universe model 
(UM) as a flat expanding FLRW one. 

The ratio Nf^ (fx/A (y)) /Nfv (fx) in Eq.E]is known 
as the amplification bias (lTurnerj|l980l . [Turner et al.l Il984l . 
IFukugita fc Turnerj|l99lh . It is introduced to take into ac¬ 
count a favourable bias when estimating the lensing proba¬ 
bility in a flux-limited sample induced by the amplification 
phenomenon. The amplification may lead to the inclusion of 
sources in a flux-limited sample that are intrinsically fainter 
than the flux limit but have undergone a gravitational lens¬ 
ing amplification. Since the ratio of the DNCF in Eq. |6]at 
two different X-ray flux levels (i.e. fx and fx/A (y)) is likely 
to be slightly dependent on the source redshift, it would cer¬ 
tainly be more accurate to consider the redshift dependence 
of the DNCF. However, at the time of this work, the XXL 
redshift are not yet available and we do not have access to 
a source sample large enough to constrain the redshift de¬ 
pendence of the DNCF (the source sample used is described 
in Section o. Consequently, we consider the value of the 
DNCF ratios averaged over the whole redshift range. Be¬ 
cause of the presence of A (y) in the amplification bias, the 
expression of T, {za, Zd, fx,a) is intrinsically linked to the 
amplification map of the deflector defining the amplification 
as a function of the source position, itself determined by 
the mass distribution of the deflector. It is out of the scope 
of this paper to explain this dependence in deta il and the 
reader may refer to e.g. iHezaveh fc Hold^ (I 2 OIII ') for some 
further description. Nevertheless, because of the presence of 
the total amplification A in its argument, the calculation of 
the amplification bias in Eq.[S] necessitates the knowledge of 
Nfj^ (fx) for sources fainter than the survey limiting flux. 
This may be estimated either by extrapolating ifx) 
or by considering a parent source sample accessing fainter 
fluxes. In this work, we will consider a deeper source sample 
described in Section [3] 

Finally, Sy in Eq. [^represents the area in the y plane, 
centred on the deflector, in which the presence of a source 
leads to a lensing event, i.e. the lensing cross section, nor¬ 
malised to the scale factor. 

The lensing event may be defined in different ways. It 
may refer to the formation of multiple images or to the for- 
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mation of a given number of lensed images, or the forma¬ 
tion of lensed images with an angular separation sufficiently 
large to be detected in a survey. Depending on the defini¬ 
tion adopted for the lensing event, the integration domain 
Sy differs and, consequently, this leads to a different defini¬ 
tion of the lensing cross section E in Eq. [6]and of the lensing 
optical depth in Eq. [S] In most cases, the integration in Eq. 
[6] must be performed numerically. 

When modelling the deflector mass distribution by 
means of an SIS profile, the lensing events can lead to the for¬ 
mation of two images at maximum. We define tsis {zs, fx) 
as the probability for a source to be multiply imaged when 
the deflectors are modelled by an SIS mass distribution. 
TSIS {zs, fx) is thus calculated by inserting Eqs. ID [S] [6] 
and [7] into Eq. [3] and by performing the integration over 
(7. The integration over a can be performed analytica lly, 
leading to the follo wing expression ilTurner et al.|[l9^ . or 
iMitchell et al.l[2005l for derivation using the VDF) 


TSIS (Zg, /x) = 'h 


/■ 

^0 


r((« + 4)//i) 
^ r(a//I) 


cdt 


(1 -I- Zd) - - EsiS (Zs,Zd, fx,(T*) dZd- 

dzd 


( 8 ) 


Esis (zg, Zd, fx, <T*) represents the lensing cross section de¬ 
fined by Eq. [6] for a — a*, when considering the area Sy for 
the case of an SIS deflector. The integration in Eq. |8] must 
be performed numerically. 

The introduction of an internal ellipticity in the mass 
distribution used to model the deflectors, allows to account 
for the formation of more than two lensed images. This is 
the reason why the Singular Isother mal Ellipsoid (SIE) mas s 
distribution has been introduced (|Kormann et al.l Il994h . 
The ellipticity parameter in the SIE profile is the axis ratio 
q of the deflector mass projected on the deflector plane and 
the mass distribution parameters M' are now a and q. As 
for the case of the SIS deflector, we can define the probabil¬ 
ity tsie {zg,fx) for a source to be lensed with the forma¬ 
tion of multiple images, irrespectively of the number of the 
lensed images, when mo delling the deflector p opulation with 
SIE mass distributions dHuterer et ahluOPSh . On the other 
hand, we may also define the probability tsie.i {zg, fx) for 
a source to be lensed with the formation of i images (with 
i = 2, 3 or 4). As the lensing cross section now depends on 
the axis ratio q of the deflector (through the dependence of 
Sy in Eq.|6}, in order to calculate tsie {zg, fx) we must also 
integrate Eq. [T] over the axis ratio q, using an appropriate 
probability distribution. Furthermore, the deflector density 
function in Eqs. |3] and |4] must take into account how the de¬ 
flector population is distributed as a function of both a and 
q, and thus we have to introduce the dependence of ud on 
the axis ratio q. In Eq. [4] the density no (cr, q) of deflectors 
with a central velocity dispersion in the range [cr, a -\- da\ 
and an axis ratio in the range [q, q -I- dq[ may be expressed 
as 


no (cr, q) = (cr) dq\^ {a, q) da 
= (a) dq (q) dadq, 


( 9 ) 


where dq^^r (cr, q) represents the normalised axis ratio distri¬ 
bution for the deflectors with a central velocity dispersion 
a, and dq (q) is the marginal normalised distribution as a 
function of q. The last equality arises if we assume that the 
deflector distributions as a function of a and q are mutu¬ 


ally independent 'ho- (cr) is given by Eq. [S] as in the case 
of the SIS mass model. There are strong evidences from the 
study of various gravitational lens samples, that elliptical 
galaxy isophotes and the mass distributio n ellipticities are 
align ed a nd have well correl ated values fsee lKoopmans et al.l 
I2OO6I and ISluse et al.l |2012| for independent confirmations). 
Furthermore, as previously mentioned, there is no evidence 
for strong evolution effects in the VDF of early-type galax¬ 
ies from lensing surveys. Thanks to these two observational 
facts, dq (q) can be estimated from the distribution of the 
isophotes of early-type galaxies as measured in the local uni- 
ve rse. We therefore u se the axis ratio distribution measured 
bv IChoi et al.l (l2007l l from a sample of elliptical galaxies in 
the local universe. 

TSIE {zg, fx) can thus be calculated by inserting Eqs. 
[ 5 ] and |9] into Eq. [T] and by integrating over cr and q. Using 
Eqs.[6]and[7l and performing the analytical integration over 
cr, TSIE can thus be expressed as 


TSIE {Zg, fx) = *!>* 


ri{a + 4)/p) 
r{a/p) 


Zs 1 

Jdq (q) Esie {Zg,Zd,fx,at,q) dqdzd, 


( 10 ) 


where Esie (zg, Zd, fx, at, q) represents the lensing cross 
section calculated through numerical integration of Eq. [b] 
when considering an SIE deflector with an axis ratio q and 
a central velocity dispersion cr*. Similarly, the probability 
TsiE,i {zg, fx) for a source to be lensed with the formation 
of i images, when modelling the deflectors with SIE profiles, 
can be calculated using EsiE,i in the previous relation and 
considering in Sy only the area in which a source should be 
located in order to lead to the formation of i lensed images. 

We have developed Matlab toolboxes and libraries to 
perform the numerical integration in the expressions of 
TSIS {zg, fx) and tsie (zg, fx) from Eqs. |8] and [TO] as well 
as for the calculation of TsiE,i {zg, fx), taking into account 
that lensed images angularly too close to each other can¬ 
not be resolved in the survey. The numerical integrations 
are made in two steps. First, we create a data base of the 
lensing cross sections in Eq. |6] considering feo = 1, ranging 
over all possible values of fx and over the ratio QmigjQE, 
where 9mig represents the smallest angular separation for 
which point-like images of same brightness are resolved in 
the survey, and also over q for the SIE case. The integration 
of the double integral in Eq.[3]is performed by plain Monte- 
Carlo integration, where we randomly generate ~ 10® source 
positions y, calculate the position and amplification of the 
lensed images, which contribute to the integral only if the 
multiple images can be resolved and brighter than the X-ray 
limiting flux of the survey. In the second step, the lensing op¬ 
tical depths are calculated by integrating Eqs. [H] and [TU] using 
trapezoidal integration and the database of lensing cross sec¬ 
tions. We have thus all the tools needed for the calculation 
of the lensing optical depth of a source with known redshift 


^ Although not strictly justified, this assumption is made because 
of the lack of observational constraints for the distribution of 
early-type galaxies in the (tr, q) plane, as well as for calculation 
time consideration. 
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and apparent X-ray flux, when considering a population of 
deflectors modelled with SIS or SIE mass distributions. 

Similar expressions for the lensing optical depths such 
as Tsis {zs, fx), tsie {zs, fx) and tsiea (zs, fx) (Eqs.[8]and 
1101 respectively), as well as their differential contribution 
as a function of the deflector redshift dT jdzd, have already 
been derived and used for the analysis of statistical samples 


fe.e. [Turner et al.lll984 iTurneilllQQCl. iFukuffita et al.l 

199C, 

Fukugita & Turnei[ll991. Kochanekll992l.lSurdei et al. 

1992, 

Maoz & Rixlll994 Cen et al.l 

1994. Keetonll2002l. Chae et al. 

200llChaell200.3l.b£ruri et al 

20121). or to studv the DODula- 

tion of deflectors 

Keeton et al. 19981. Keeton & Kochanek 

19981. 

Keeton 1998 

. Kochanek et al.l 2000l. Ofek et al.l 20031. 

Chad 

2 OO 3 I. IChaell2010lj. some of the work having consid- 


ered the ellipticity of the mass distrib ution as well as the 
effect of external sheer on the statisti cs dHuterer et ahllioosl . 
lOeuri fc Marshaiill2010l . [Oguri et al.ll2012l b 

The various expressions derived for the optical depth r 
and its differential contribution dr/dzd are considered for a 
single source. This may be averaged over the whole detected 
population to derive the mean optical depth through the 
sample as well as the expected redshift distributions of the 
lensed sources and deflectors. To average over the popula¬ 
tion of sources, some previous works have made use of the 
source luminosity function, and integrate over the absolute 
magnitude of the sourc e population (e.g. lOeuri fc Marshall 
I 2 OI 0 I . [Oeuri et al]l2012l i which, for the integration process, 
necessitates the choice of a universe model. 

In the next subsection, we propose a slightly different 
formulation that allows to average any function over the en¬ 
tire population of sources detected in a survey, where the 
averaging is done based on the observables, i.e. the distribu¬ 
tion in the {zs,fx) plane of the detected sources. Because 
the formalism is based directly on the observed distribution 
of sources, it naturally takes into account the detection bi¬ 
ases of the sources. We apply this method to derive useful 
expressions such as the average lensing optical depth in a 
sample, the redshift distributions of the lensed sources and 
of the deflectors effectively leading to the formation of mul¬ 
tiple lensed images. 


2.2 Joint source probability density and fraction 
of multiply-imaged sources 

Let us consider a survey characterised by its limiting flux, 
different biases in the source detection procedure and its 
angular coverage. Each detected source is characterised by 
its redshift Zs and its flux fx in the selected spectral band 
(in the present case, the [0.5 — 2] keV band). The detection 
of a source within the survey may be considered as a ran¬ 
dom event with respect to the continuous random variables 
associated with the source redshift and the X-ray flux, re¬ 
spectively. 

We can define a probability P{za,fx) that a source 
detected in the survey is characterised by an observed red¬ 
shift and flux in the range [za, Zs + dza] and [fx, fx + dfx], 
respectively. We may then define the joint probability den¬ 
sity dobs (zs, fx) spanning over the {za, fx) plane, associated 
with this random event. P {za, fx) and dobs {za, fx) are re¬ 


lated through 

P (Za, fx) = dobs {Za, fx) dZadfx- (H) 

The random variables Za and fx associated with a de¬ 
tected source follow the joint distribution described by 
doba{za,fx)- This functiou contains all the information 
about the survey and implicitly takes into account the de¬ 
tection biases. For a sufficiently large number of detected 
sources, the joint probability density doba Gs,fx) may be 
directly estimated from the detected source population, by 
calculating a smoothed histogram of the source distribu¬ 
tion in the {za, fx) plane, normalised by the total num¬ 
ber Nagn of sources detected within the survey. Ideally, 
we would have liked to define the source joint probability 
density doba {za,fx,r) in the {za,fx,r) space, where r rep¬ 
resents the SDSS r-magnitude of the optical counterpart (s) 
of the X-ray sources. However, because of the small num¬ 
ber of sources observed in our reference sample (see Section 
0. the quantity doba {za,fx,r) could hardly be accurately 
determined. When analysing the XXL sample however, the 
number of detected sources should be sufficient to constrain 
doba in the 3D-space. We shall therefore postpone such a 
more detailed study until all optical counterparts of the XXL 
X-ray sources have been identified. 

The normalised marginal probability density distribu¬ 
tions associated with the random variables Za and fx are 
closely related to the observations. Indeed, the marginal den¬ 
sity distribution obtained by integating doba{za,fx) over 
Za or fx represent the normalised source distribution as 
a function of the flux ifx) /Nagn and the redshift 
Nza [za) /Nagn, respectively. 

The use of doba{za,fx) as a weighing function when 
performing the integration over the entire population of de¬ 
tected sources allows to calculate the expected mean value of 
any function of the random variables Zs and fx ■ The mathe¬ 
matical expectation (r) of the lensing optical depth, i.e. the 
fraction of sources gravitationally lensed within the detected 
population, can be calculated by integating r(za,fx) over 
the (za,fx) plane, weighing with doba {za,fx), i.e. 

(t) = JJ T{Za,fx)doba{Za,fx)dZadfx- (12) 

Using the expression of Tsis{za,fx), TsiE{za,fx) or 
rsiE,i (za, fx) given by Eqs. |8] and [10] for the calculation 
of the optical depth T{za,fx) in Eq. 1121 we are able to 
calculate the expected fraction of multiply-imaged sources, 
considering a population of deflectors modelled with either 
the SIS or SIE mass profile. In the latter case, we can also 
calculate the expected fraction of lensed sources as a func¬ 
tion of_thejiurnbCTjof lensed images. 

lOeuri et al.l 1 I 2 OO 8 I I have also derived an expression for 
the number of multiply-imaged sources using a binning o f 
the redshift-magnitude space (Eq. 12 in lOguri et al.l[2008l L 
and with the number of sources in the bins as a weighing 
factor. Their expression corresponds to the discrete equiva¬ 
lent of Eq. integrating in the redshift-magnitude space, 
rather than the {za, fx) plane. The essential difference with 
the approach used in these previous works and ours is the 
use of the source distribution in the {fx, Za) rather than us¬ 
ing the source luminosity function as a weighing factor in 
the absolute magnitude-redshift space. Our approach allows 
to account directly for the detection bias of the sources and 
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does not necessitate any assumption of a universe model for 
the calculation of the weighing factor (which is necessary 
when using absolute magnitudes and the source luminosity 
function). 

2.3 Normalised redshift distributions 

It is now straightforward to establish the normalised distri¬ 
bution wz^ (zd) of the deflector redshifts expected in the 
XXL field 



(13) 

where the differential contribution dr/dzd of the redshift Zd 
to a source leasing optical depth is given by Eq. [T] 

Similarly, the normalised redshift distribution wz^ (zs) 
of the leased sources is given by 

Wz, (Zs) = Jt {Zs,fx) dobs {Zs, fx) dfx- (14) 

lOeuri et alj ll2012li have derived by different means an 
expression for wz^ (zd) an d applied it to the sample of the 
SDSSQLS (see Eq. 23 in lOguri et al.l[2012l i. This estima¬ 
tion was done by binning the redshift-magnitude plane and 
using the number of sources in the different bins as a weigh¬ 
ing factor. Their relation corresponds to the discrete equiv¬ 
alent of Eq. m where the integration runs over the redshift- 
magnit ude plane, rather t han t he {zs,fx) one. 

In lOguri fc Marshalll (l2010l l. the authors have derived a 
different expressi on for the redshift distrib ution of the leased 
sources (Eq. 7 in lOguri fc Marshallll201^ 1. equivalent to Eq. 
M except that the integration runs over the absolute mag¬ 
nitude and the weighing function used in the integration 
corresponds to the expression of the joint probability den¬ 
sity expressed in terms of t he source luminosity f unction. 

Finally, let us note that iMitchell et aP (l2005l l have also 
derived an estimation for the redshift distribution of the 
leased sources in the CLASS survey, assuming the DNCF 
of the sources to be expressed as a single power-law expres¬ 
sion, i.e. the amplification bias being thus constant for each 
source. However, this assumption is very restrictive as QSOs 
usually show a DNCF that presents a break at a critical 
magnitude, and needs to be modelled by a double power 
law expression. 

The use of dobsizsjfx) allows an easy calculation of 
(r) through the survey, and the calculation of the expected 
normalised distributions wz^ {Zs) and wz^ iz:d), as a func¬ 
tion of the redshift of the lensed sources and of the deflec¬ 
tors, respectively. These distributions are calculated without 
any assumption about the source population, as dobs {zs, fx) 
may be directly estimated from the observed data, naturally 
including the observational biases. 

In the next Section, we present the observational con¬ 
straints used for the estimation of the joint probability den¬ 
sity of the XXL survey. 

3 OBSERVATIONAL CONSTRAINTS AND 

Dobs DETERMINATION 

Ideally, the joint probability density of the XXL sources 
should be derived by doing a histogram of the XXL 


sources in the redshift-X-ray flux plane. However, the cross¬ 
correlation between the X-ray and (near-)optical data and 
the identification of the source type have not yet been per¬ 
formed. The XXL X-ray sources and their associated multi¬ 
wavelength data have characteristics similar to those of one 
of its sub-fields, the XMM-Large Scale Structure (XMM- 
LSS), covering 10.9 out of the 44.2 square degrees effec¬ 
tively covered by the XXL survey. The XMM-LSS X-ray 
sources and their associa t ed m ulti-wavelength data, pre¬ 
sented in IChiappetti etld] ll2013tl . have optical counterparts 
taken from the CFHTLS W1 catalogue, down to the lim¬ 
iting magnitudes i' ~ 25, r' ~ 25, before correction for 
Galactic extinction 0. The X-ray source classification and 
propert ies of the XMM-LS S field have already been deter¬ 
mined llMelnvk et ahll^OlSl j. The survey limiting flux in the 
soft band is F[Q,b- 2 ]keV — 3 x 10“^® erg s~^ cm~^ (with a 
detection probability of 0.5, as defined in lElviv et al .11201^ 1. 
Nevertheless, at the time of this work, the redshift estimate 
in the XMM-LSS sample is still ongoing. 

Therefore, to retrieve the properties of the expected 
population to be detected within the XXL survey and their 
optical counterparts, we used data from a deeper held (~ 60 
ks exposure, compared to ~ 10 ks for XMM-LSS): the 
XMM-COSMOS held, to which we apply the hux cuts of 
the XXL in the X-ray and the optical, that we assume to 
be similar to those of the XMM-LSS. The XMM-COSMOS 
survey, cover ing a contig u ous h eld of 2 square degrees is 
described in iBrusa et al.l (l2010l ). and has a limiting hux 
of E[o,5_2]fcev — 5 X 10“^® erg s“^ cm~^ (~ 1.4 x 10“^® 
considering the hu x with with 50% effe ctive area cover¬ 
age, see Fig. 6 in ICappelluti et ^ l2007t) . The hux limit 
at 50% effective area coverage for the XMM-LSS and the 
COSMOS are in good agreement with the exposure time 
ratios as 3 x 10~^®/1.4 x 10“^® ~ 60fcs/10fcs. Their op¬ 
tical catalogue with which the cross correlation was per¬ 
formed contains sources detected in at least one of the Sub¬ 
aru bands (b, v, g, r, i, z) down to an AB magnitude limit 
of ~ 27. The COSMOS sample is almost complete down 
to F^o. 5 - 2 ]keV — 3 X 10~^ ® erg s~^ cm~^ (X-ray detection 
probability of 0.98, cfr ICappelluti et (l2007l j. Fig. 6) 
and 98% of the X-ray sources have an optical counter part. 

The joint probability density of the XXL survey is ex¬ 
pected to be quite similar to that of the XMM-COSMOS 
held if we apply to this survey the same X-ray and optical 
hux cuts and take into account the different probability of 
detection for sources fainter than 2 x 10~^"^ erg cm“^ s“^. 
Consequently, we hrst estimate the COSMOS joint probabil¬ 
ity density dcosMOS {zs, fx) from a smoothed histogram of 
the sources from the XMM-COSMOS in the {zs, fx) plane, 
where we apply the X-ray and r-band cutohs of the XXL. 
We then take into account the XXL detection probability to 
determine dobs (zsjx). 

In the left panel of Fig. [T] we have represented the 
XMM-COSMOS sources in the (zs^fx) plane, with fluxes 
presented along a logarithmic scale. We have only considered 

^ After correction, the limiting magnitudes of the CFHTLS 
W1 catalogue are (expressed in AB magnitudes and considering 
a 5 sigma signal-to-noise ratio) r ~ 24.8 and i ~ 24.5. See 
http: //xxlmultiwave.pbworks.com/w/page/54613008/0ptical 
and http://terapix.iap.fr/rubrique.php?id_rubrique=268 for 
details. 
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Figure 1. The left panel displays the distribution in the (zs, fx) plane of the XMM-COSMOS source population llBrusa et al.|[2oT3l , 
restricted to F^o.5—2]keV > 3 x 10“^^ erg cm~^ s“^ in the [0.5 — 2] keV band and with an optical counterpart brighter than r = 25, for 
which a reds hift estimate is ava ilable. This X-ray band flux limit corresponds to a detection probability of 0.5 in the XMM-LSS fleld, 
as defined in lElviv et al.l ll2012l') . The X-ray flux is shown on a logarithmic scale. The right panel displ ays the joint prob ability density 
function (zs^fx) in the {zs^fx) plane, corresponding to the XMM-COSMOS source population ijBrusa et al.l[2Ql^ . restricted to 
the X-ray band flux limit of the XMM-LSS fields. 


sources with a flux larger than -F[o. 5 - 2 ]fceV — 3 X 10”^® erg 
cm~^ in the [0.5 — 2] keV band and with an optical coun¬ 
terpart brighter than r = 25 in the SDSS r-band. We have 
rejected all X-ray sources for which there was no redshift 
estimate. Whenever available, spectroscopic redshifts were 
preferred to photometric ones. The XMM-COSMOS cata¬ 
logue presented in iBrusa et al.l (l2010l l contains 1797 sources 
in the 2 square degrees. After applying the X-ray and r-band 
cut-off (and excluding sources with no available r magni¬ 
tude), there are 630 sources of which 6 are excluded because 
of no available redshift. The final restricted COSMOS source 
sample contains 624 sources. The density of sources in the 
{zs,fx) plane is larger for the fainter fluxes, near redshift 
2 ~ 1. At any given redshift, there are more sources with 
a lower flux. Finally, let us also point out the absence of 
sources with a high flux at high redshift. 

In order to estimate the joint probability density 
dcosMOS {zb, fx) of the COSMOS sources, we have calcu¬ 
lated a smoothed histogram of the COSMOS source dis¬ 
tribution with the XMM cut of in the X-ray and r-band 
displayed on the left panel of Fig. [T] For the convenience 
of the developed software, dcosMOS [zs, fx) has been de¬ 
rived using a logarithmic scale for the X-ray fluxes. We have 
considered redshift intervals of 0.375 with bin centers sepa¬ 
rated by 0.0625 and logarithmic magnitude intervals of 0.3 
separated by bins of 0.05. The derived COSMOS joint prob¬ 
ability density is shown in the right panel of Fig. [T] We have 
intentionally kept the same axis as in the left panel in order 
to clearly identify the similarities between the two figures. 
The grey scale indicates the values of dcosMOS [zs, fx)', 
the darker the grey, the higher the probability of finding a 
source. The normalisation factor of dcosMOS (za, fx) is the 
number of sources (624) detected in the 2 square degrees of 
the COSMOS field, restricted to the XMM-LSS cutoffs. 

To take into account the detection probability of the 
XXL survey, we multiply the dcosMOS (za, fx) in the right 


panel of Fig. □ by the XMM-LS S detection proba bility as a 
function of the flux taken from lElviv et ah] (l2012l l. Fig. 10, 
the obtained distribution being doba {za,fx)- The normali¬ 
sation factor of this distribution represents the number of 
sources (592.5) with Fio. 5 - 2 ]keV 3 x 10“^® erg cm~^ s“^ 
and r < 25 that would be detected within the XXL sur¬ 
vey in the 2 sq. deg. of the XMM-COSMOS field. From this 
normalisation factor, we may thus estimate the number of 
sources to be detected in the XMM-LSS and XXL-field with 
similar characteristics. The expected numbers of sources in 
the different surveys are summarised in Table [T] Assuming 
the final XXL catalogue to have similar characteristics as 
those of the XMM-LSS, this estimate of doba {za,fx) is as¬ 
sumed to be valid for both the XMM-LSS and the XXL 
fields. 

As a reliability test of the derived joint probability den¬ 
sity doba {za, fx), we Verify its ability to represent the prop¬ 
erties of the observed population of AGNs. In the left panel 
of Fig. [21 we have represented as a continuous grey line the 
observed cumulative distribution as a function of the redshift 
of the XMM-COSMOS sources with an X-ray flux larger 
than F’[o.5-2]fceV = 3 x 10“^® erg cm“^ s“^ in the [0.5 — 2[ 
keV band, with r < 25 and a redshift estimation. 

In the same figure, we have represented with a dashed 
dark-grey line the cumulative source redshift distribution de¬ 
rived from the COSMOS joint probability distribution when 
not considering the detection probability of the XXL survey. 
The cumulative redshift distribution is obtained by integrat¬ 
ing dcosMOS (za, fx) over fx and by multiplying the num¬ 
ber of AGNs detected per square degree in the COSMOS 
field, where we have applied the XXL cutoffs. 

The cumulative redshift distribution inferred from the 
distribution dcosMOS [za, fx) reproduces very well the ob¬ 
served cumulative source redshift distribution. 

On the same figure, we have represented with a dark 
continuous line the cumulative redshift distribution derived 
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Figure 2. Left panel : Cumulative function of the number density N {zg) of the sources as a function of the redshift for the XMM- 
COSMOS sources with F[o.5—2]keV ^ 3 X 10“^® erg cm“^ s“^ and r < 25. We have also represented the cumulative redshift distribution 
inferred from the marginal distribution of the joint probability density function dcoSMOS fx) and d^bs {zs, fx), with and without 
considering the XXL detection probability detection (dashed and continuous black curves, respectively). Right panel: Differential number 
counts of the XMM-COSMOS sources as a function of the flux in the soft X-ray band, tor all the COSMOS sources, which are used for 
the amplification bias calculation. We have only represented the X-ray flux range accessible to the XXL survey. We have also represented 
the marginal distribution obtained from the joint probability density function dcoSMOS {^s,fx)^ as well as the DNCF fit used for the 
calculation of the amplification bias. 


from the final dobs {zs, fx) when considering the detection 
probability of the XMM-LSS survey. We remark a very good 
agreement at low redshift between this distribution and the 
two previous ones, and a significant deviation above redshift 
~ 1 due to the fact that the XXL survey does not detect all 
of the fainter sources, more numerous at higher redshifts. 

In the right panel of Fig. [J] we have represented the 
differential number counts function (DNCF) of the XMM- 
COSMOS sources as a function of their X-ray soft band 
flux. The flux is shown along a logarithmic scale and the 
error bars are estimated considering a Poisson noise. We 
have here considered all the COSMOS sources brighter than 
F[Q b- 2 ]keV = 1 X 10“^® erg cm“^ s~^. This DNCF is used to 
estimate Nf^ ifx) needed when calculating the amplifica¬ 
tion bias in Eq. |6] As previously mentioned, the calculation 
of the lensing optical depth necessitates the knowledge of the 
DNCF for fluxes fainter than the survey detection limit. We 
have therefore fitted the observed XMM-COSMOS DNCF 
logio ifx)) with a 3’’'^ order polynomial in the flux range 
covered by the XMM-COSMOS and extrapolated the data 
linearly outside this flux range. The fit is shown as a contin¬ 
uous dark-grey line in the right panel of Fig [2] This poly¬ 
nomial form of the DNCF is used for the calculation of the 
amplification bias in our simulations. 

On the same figure, we have represented with a dashed 
line the DNCF per unit of solid angle derived from the joint 
probability density function dobs (zs, fx), where we have re¬ 
stricted the COSMOS sources to the XXL cutoffs and taken 
into account the detection probability of the survey. This is 
obtained by integrating dobs [zs, fx) over Zs and by multi¬ 
plying by the number Nagn/^ of AGNs per square degree. 
This curve is matching fairly well the fitted curve for the 
brighter sources and is lower for the fainter sources, where 


the XXL survey only detects part of the sources and where 
the fraction of excluded sources with r > 25 is higher. 

The small size of the source sample in the XMM- 
COSMOS field at low redshift (due to the small survey 
angular coverage) does not permit a better estimate of 
dobs (zs, fx) at these low redshift values, because of a large 
scatter in the observed data at low redshift, especially for 
the brighter sources. Nevertheless, as these very bright and 
low redshift sources are very rare in the survey and as their 
lensing probability is very small (because of their very low 
redshift), this does not have a large impact on our simula¬ 
tions. These problems will be reduced in the XXL survey, for 
which the angular coverage is ~ 20 times larger than that 
of the XMM-COSMOS field. This confirms, as mentioned 
earlier, that the COSMOS sample is not large enough to de¬ 
termine the distribution dobs {zs, fx,r) of the sources in the 
3-D space (zs,fx,r). 


4 RESULTS 

4.1 Mean lensing optical depth 

Using the joint probability density dobs (zs, fx) described in 
Section[3l we have computed the mean lensing optical depth 
(r) for the XXL survey, integrating numerically Eq. 1121 As 
dobs {zs, fx) is assumed to be identical for the XXL and the 
XMM-LSS surveys, the results are valid for both surveys. 

To calculate the lensing optical depth, we modelled the 
deflector population with SIS and SIE mass distributions 
(Eq. [8]and Eg. 1101) . We have computed the average optical 
depth for different values of the minimum image separa¬ 
tion 6mis resolvable at optical wavelengths. When calculat¬ 
ing the cross section, we have considered the detection of 
the lensed images to be achievable down to an angular sep- 
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Figure 3. Average lensing probability ( t ) calculated for the XXL 
and XMM-LSS fields as a function of the minimum angular sepa¬ 
ration 6mis resolvable by the survey in the optical domain, when 
modelling the deflector population with the SIS and SIE mass dis¬ 
tributions. For comparison, we have also computed the case of the 
late-type galaxies, modelled by means of SIS mass distributions. 


aration 9mis (independently of their relative amplification). 
Although this is not strictly accurate, we have made this 
assumption for the following reason. The regions contribut¬ 
ing the most to the lensing cross section are those where 
the source is located close to the caustic curves, as these 
are the most amplified and thus beneht the most from the 
amplification bias. For these configurations, the lensed im¬ 
ages that are the closest to each other (and also the bright¬ 
est) are those formed on each side of the tangential critical 
curv e. These lensed image s have a very similar amplifica¬ 
tion llKormann et al.lll994h . Consequently, for the detection 
of these lensed images, the critical parameter is the smallest 
angular distance under which these point-like images cannot 
be resolved, independently of their brightness. In practice 
however, for lensed images close to each other, the mini¬ 
mum angular distance at which the point-like images can 
be disentangled is dependent on the flux difference between 
the images, especially at very small angular distances. For 
the analysis of the final XXL sample, the angular selection 
function (characterizing the smallest angular distance de¬ 
tectable as a function of the relative amplification of the 
lensed images) will have to be determined precisely. 

Fig. |3] displays the behaviour of the average optical 
depth (r) as a function of the minimal image separation 
dmis resolvable in the optical survey. We display the results 
when modelling the deflector population with either SIE or 
SIS mass profiles. In the former case, we have calculated the 
total probability of having a multiply-imaged source, inde¬ 
pendently of the number of the lensed images. 

For a perfect instrument (i.e. Omis = 0”), when consid¬ 
ering deflectors modelled with SIE mass distributions, we 
find an average lensing probability (tsie) = 1.698 x 10“®. 
When modelling the deflectors with an SIS mass distribu¬ 
tion, the mean lensing probability is (rsis) ~ 1.788 x 10“®. 
Modelling the deflectors with SIS mass distributions thus 
leads to a slightly larger average lensing probability by 
~ 5%. The reason for this is that the SIE mass distribu¬ 


tion was introduced in order to preserve the projected mass 
inside a same area but does not conserve the geometrical 
cross section nor the amplification probability distribution. 
In other words, the area inside the caustic curve of an SIS 
deflector is always larger or equal to that included inside 
the caustic curves of an SIE deflector. When averaging the 
SIE lensing cross section over the deflector ellipticity distri¬ 
bution, this leads to an equivalent SIE lensing cross section 
smaller than that of the SIS case. Furthermore, for these two 
types of deflectors, the probability of producing a multiply- 
imaged source with a given total amplification slightly dif¬ 
fers. The impact of the amplihcation bias is thus different for 
the two deflector models. Consequently, the slight differences 
between (tsis) and (tsie) depend on the DNCF as a func¬ 
tion of fx for the source population, which varies from one 
survey to another. The overestimate of (tsis) thus has to be 
estimated independently for each survey. This boost of the 
av erage lensing o p tical depth in the SIS model was studied 
by iHuterer et al.l ll2005h who concluded that the ellipticity 
in the deflector mass distribution decreases the mean lensing 
optical depth. A boost of the SIS model by a few percent 
may be expected in a survey with sources showing a steep 
DNCF. 

Both (tsie) and (rsis) decrease for increasing values 
of the parameter 9mis- When considering a finite resolution 
of the instrument, some of the lensed images formed are 
angularly too close to each other and are detected as a single 
point-like object. Consequently, the probability of detecting 
the lensed images decreases as 9mis increases. 

In the present case, the multiple lensed images will be 
searched for among the optical counterparts of the point¬ 
like X-ray sources, because of the better angular resolution 
in the optical domain. The ground-based observations are 
limited by the atmospheric seeing. In the northern XXL- 
helds , the CFHT in the r-band have a typical seeing of 
0.7” llSalmon et al.l 1^00911 and for the southern hemisphere 
the typical seeing with the Blanco telescope is ~ 0.9” 
JPesai et al.l I 2 OI 2 I 1 . We thus considered our full width at 
half maximum of the PSF to be homogeneous over the en¬ 
tire sample and equal to ~ 0.9”. 

Thanks to point spread function fitting techniques, we 
can hope to resolve multiple point-like lensed images down 
to half the full width at half maximum for lensed images 
with a same amplification, which constitute the conhgura- 
tions contributing the most to the lensing cross section as 
previously explained. Consequently, the typical 9mis value 
achievable is expected to be 9mis ~ 0.45” (in practice, as 
stressed previously, 9mis depends of the relative amplifica¬ 
tion of the lensed images, which will have to be taken into ac¬ 
count when analyzing the final XXL sample). For this value 
we have (tsie) = 1.384 x 10“® and (tsis) = 1.489 x 10“^. 
The slight overestimate of the SIS mean lensing value rela¬ 
tively to that of the SIE model thus increases with 9mis and 
reaches ~ 10% for 9mis = 0.45”. 

For comparison, we have computed the evolution of the 
mean optical depth as a function of 9mis for the population 
of late-type galaxies modelled by means of an SIS mass dis¬ 
tribution. As for the case of the early-type galaxies, we have 
considered the comoving density of late-type galaxies to be 
constant with the redsh ift and we h ave used the VDF pa¬ 
rameters determined by IChael (l2010h in the local universe. 
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Table 1. Estimate of (1) the number of sources with an X-ray flux greater than 3 x erg cm“^ s“^ in the [0.5 — 2]keV range and 

r < 25 and (2) the number of detected multiply-imaged sources in the XMM-COSMOS, XMM-LSS and XXL surveys. The estimates 
of the number of sources are extrapolated from the XMM-COSMOS catalogue taking into account the angular coverage of the different 
surveys. For the estimate of the number of multiply-imaged sources, the numbers in parentheses correspond to the SIS case. 


Survey 

Coverage 

(deg2) 

Number of sources 

Lensed sources 
{dmis = 0”) 

SIE (SIS) 

Lensed sources 
{dmis=0A5”) 
SIE (SIS) 

>2 images 

{emis= 0 ”) 

>2 images 
{dmis = 0.45”) 

XMM-COSMOS 

2 

621 

- 

- 

- 

- 

XMM-COSMOS 

(with XXL Pdetection) 

2 

529.5 

- 

- 

- 

- 

XMM-LSS 

10.9 

2885 

5 (5) 

4(4) 

0 

0 

XXL 

44.2 

11701 

20 (21) 

16 (17) 

1 

1 




{t} 

{Bmis = 0”) 

{t} 

{Bmis =0.45”) 

(r> 2 }/(r) 

{dmis =0”) 

(t>2) / (t) 
{Omis = 0.45”) 

SIE 



1.698 X 10“® 

1.384 X 10“3 

0.0718 

0.0564 

SIS 



1.788 X 10-3 

1.489 X 10-3 

- 

- 

SIS (Late) 



5.719 X lO-'^ 

1.629 X 10-3 

- 

- 


i.e. 

[$*,o-*,a,/3] = [66 X 10"® h®Mpc■^ 91.5 kms"\0.69, 2.10] . 

The results are displayed in Fig.O When considering a per¬ 
fect instrument {9mis = 0”), the average lensing optical 
depth associated to the late-type galaxies is about a third 
of that of the early-type ones (i.e. (r) = 5.719 x 10"'^). 
Nevertheless, the decrease of (r) with Omis is steeper than 
in the early-type galaxy case. Indeed, late-type galaxies are 
less massive and lead to smaller typical angular separations 
of the multiple lensed images, which are not disentangled 
in the seeing limited images. For 6mia = 0.45”, the average 
lensing optical depth due to late-type galaxies is found to 
be (r) = 1.629 x 10“®, three orders of magnitude lower than 
that of the early-type galaxies. The expected contribution of 
the late-type galaxies in our sample of gravitationally lensed 
sources is thus negligible, which validates our assumption of 
only considering the population of early-type galaxies as the 
deflectors for the XXL lensed source sample. 

Our estimation of the contribution of late-type galax¬ 
ies to the lensed sources is surprizingly low compared to 
the observed fraction of late-type lenses in existing samples. 
Furthermore, it suggests that late-type lenses with image 
separations larger than 0.5” are extremely rare, which is 
also in contradiction with observed samples. Indeed, 2 out 
of the 13 CLASS lenses from the statistical sample are likely 
to be produced by l ate-type galaxies (B0218 and B1933, see 
iBrowiie et al.l 1^0031 ) and out of the 26 lensed QSOs of the 
SDSSQLS statistical sample, one is possibly due to a late- 
type galaxy (J1313) and has an angular separation larger 
than 1”. We do not fully understand the reason for these 
discrepancies. A possible cause of error is the effect of the 
lens environment, not considered in this work, which may 


lead to an additional gravitational shear. Huterer et al. 2005 
have shown that the external shear broadens the distribu¬ 
tion of angular separation between the lensed-images with¬ 
out changing its average value. It may therefore increase 
the fraction of events with an angular separation larger that 
1”. In the CLASS and SDSSQLS statistical samples, all the 
lenses produced by late-type galaxies with and angular sep¬ 
aration larger than 1” required external shear to accurately 
model the position and relative amplification of the lensed 
images fsee lSuvu et al.l[201^ . ISluse et al.ll2012l and references 
therein). Another possible sourc e of discrep ancy is the VDF 
used for late-type galaxies (from Ic^aell2010fl which does not 
come from direct measureme nts (as in the case of the early 
type VDF. [Choi et alll2007l L It is inferred from the local 
luminosity function of late-type galaxies using the Tully- 
Fisher relationship (taking into account its dispersion) and 
assuming a conversion between the circular ve locity Vc and 
cr to be that of an SIS profile (i.e. a = u,./\/2L fchael (|201flh 
stresses that lensing statistics of late-type galaxies might 
necessitate to consider the circular velocity function (rather 
than the VDF) and consider more realistic mass distribution 
models. Furthermore, type-specific LF (and VDF) are still 
potentially biased by misclassifi c ation of the galaxy types 
and, a ccording to IPark fc Choil (|2005l '> , the work of IChael 
ll2010ll has a ~ 10% classification mismatch. We therefore 
advise to use with caution the results concerning the late- 
type population statistics. 

Table [T] summarises the number of multiply-imaged 
sources expected in the different surveys as well as the ex¬ 
pected number of detected events. Assuming the sources de¬ 
tected in the XXL survey will have the same properties as 
those in the XMM-COSMOS to which we applied the same 
flux limits in the X-ray and optical bands (except when esti- 
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Figure 4. Fraction {TsiE,i) / Psie) of the lensing events with 
a given number i = 2, 3 or 4 of images as a function of the survey 
minimum angular separation 9miat when modelling the deflector 
population with the SIE mass distribution. 3-image configurations 
are quads with 2 blended point-like images. 


mating the amplification bias) and when accounting for the 
detection probability of the XXL, we expect the detection 
of 11701 sources with i^[o. 5 - 2 ]fcey > 3 x 10“^® erg cm“^ s“^ 
and r < 25 in the 44.2 sq.deg. of the survey. Among these 
sources, we expect 21 (20 in the SIE case) to be multiply im¬ 
aged, out of which 17 (16 in the SIE case) should be detected 
assuming dmia = 0.45”. 

When calculating the lensing probability of a source to 
be lensed through Eq. 1101 considering a population of de¬ 
flectors modelled with SIE profiles, we may also calculate 
the probability rsiE,i of a source to be lensed with the for¬ 
mation of a given number i of lensed images. We have cal¬ 
culated the average value of (tsieA for the population of 
sources to be detected within the XXL and XMM-LSS fields 
using Eg. 1121 for different values of the Bmia parameter. The 
results are displayed in Fig. [4] where we have plotted as a 
function of the value of the Bmia parameter, the fraction of 
lensing events composed of i lensed images (tsieA / (tsie) 
relatively to the total average lensing probability (tsie)- 

The fraction of lensing events with the formation of two 
images is always the highest. This is a consequence o f the 
scarcity of very elliptical deflectors fsee lChoi et al.l[2007l . Fig. 
13 for the axis ratio distribution of early-type galaxies). In 
the case of a perfect instrument (i.e. Bmia =0” ), we find 
that 93% of the lensed sources are composed of 2 images and 
this fraction increases with the value of the Bmia parameter. 

For a perfect instrument, the lensed sources with forma¬ 
tion of more than 2 images is composed of quads (i.e. 4-image 
configurations). As Bmia increases, some of the 4-lensed im¬ 
age configurations, due to their too small angular separation, 
have only 3 point-like images detected (quads with 2 blended 
point-like images). Out of the 20 lensed sources formed in 
the XXL population, only one is expected to be detected 
with more than 2 images. 

The fraction of lensing systems with more than 2 images 
is rou ghly consistent with the results of lOguri fc Marshalll 
ll2010li who have calculated the expected number of gravi¬ 
tationally lensed quasars in wide-held optical surveys. For a 


Figure 5. Normalised redshift distributions of the deflectors, the 
lensed sources and all AGNs. For each distribution, we have in¬ 
dicated the median redshift value. 

survey with a limiting magnitude i = 25, these authors hnd 
a fraction of a little more than ~ 10% of quads (formation of 
4 images) against ~ 7% in our simulations. The slight differ¬ 
ence is most likely due to the different ellipticity distr ibution 
of the deflectors considered llOguri fc Marshal]||2010l consid¬ 
ers a combination of oblate and prolate three dimensional 
deflectors with a Gaussian distribution of their ellipticity) 
as well as their consideration of an additiona l exter nal shear 
due to the lens environment. iHuterer et al.l (l2005l l showed 
that the external shear increases the fraction of quads in a 
sample of lensed sources. 


4.2 Redshift distributions 

Using the joint probability density doba (za, fx) described 
in Section |3] we have computed the mathematical expec¬ 
tation of the normalised redshift distribution of the deflec¬ 
tors wz^ {zd) and of the lensed sources WZa {za), numeri¬ 
cally integrating Eos. 1131 and 1141 respectively. We have con¬ 
sidered Bmia = 0.45” and we have modelled the deflectors 
with SIS mass distributions. The normalised redshift distri¬ 
butions wzj^ (zd) and wz^ (zs) are shown as a function of 
the redshift in Fig. 0 For comparison, we have represented 
the marginal distribution as a function of the redshift of the 
joint probability density, i.e. the normalised distribution as 
a function of the redshift of all the sources (independently 
of the fact that they are being lensed). 

On this figure, we have also illustrated the median value 
of the observed redshift for the different distributions. The 
redshift distribution of the lensed sources is shifted towards 
a higher redshift compared to that for all the sources. As 
the source redshift increases, so does its geometrical lensing 
volume (i.e. the volume in which the presence of a deflector 
leads to the formation of multiple lensed images). Conse¬ 
quently, sources with a higher redshift tend to have a higher 
lensing probability. The mathematical expectation for the 
redshift moves from {zagn) — 1.11 for the entire source 
population to {za) — 1.8 for the lensed sources. 

For redshifts larger than {zagn), the distribution of the 
lensed sources does not appear as a smooth function of the 
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redshift, and we clearly see the presence of bnmps or red- 
shift ranges with a probability excess compared to a smooth 
decreasing function of the redshift. These probability over¬ 
densities of lensed source detections correspond to redshift 
ranges in which strong emission lines of the AGNs enter in 
the optical SDSS r-band, in which the optical counterparts 
are searched for; they are thus a consequence of a selec¬ 
tion bias. The presence of the emission lines in the r-band 
increases their probability of being detected, compared to 
that of a source with only a continuum-like spectrum. For 
example, the Mgll line at 279.8 nm enters the r-band Hlter in 
the redshift range 2 - 2.42, for which we see an over-density 
in the lensed source redshift distribution. 

The redshift distribution of the lensed sources corre¬ 
sponds to the probability density from which the detected 
lensed source redshift may be considered as a random event, 
thanks to the fact that, in this case all the lensed sources 
are detected, and may have their redshift estimated. Given 
a large enough sample of lensed sources, this distribution 
could be retrieved from a normalised histogram of the lensed 
sources, as a function of their redshift. 

In Fig. [5l we have also represented the normalised dis¬ 
tribution as a function of the redshift, for the deflectors 
involved in the formation of multiply-imaged sources. The 
deflector redshift median value is (zd) ~ 0.54 and the most 
probable value is below 2 : < 0.5. The contribution of deflec¬ 
tors with 2 > 1 is very small. This comforts our assump¬ 
tion of a non-evolving deflector population in the calcula¬ 
tion of the lensing optical depth, as the population involved 
is mainly located at low redshift. 

In the case of the deflector distribution, the observed 
redshift distribution of the deflectors involved in the forma¬ 
tion of multiple images of a source will be highly biased as 
mos t of them are not bright e nough to be detected. 

lOguri fc Marshall (l2010l ) have also estimated the ex¬ 
pected lens redshift distribution for lensed quasars detected 
in optical imaging surveys. The deflector redshift distribu¬ 
tions are marginally consistent, although the distribution 
fou nd in this work peaks at lower redshifts (zmax ~ 0.6 
for lOguri fc Marshalll 12^0 and Zmax ~ 0.5 in the present 
work). The expected source redshift distribution is also 
shifted towards lower redshift in our study. This difference in 
the lensed source and deflector redshift distributions comes 
from the difference in the source distribution: in the present 
work, the rather bright X-ray flux cut-off tends to reject 
sources at high redshift, that are included in the sample of 
lOguri fc Marshaiil (l201oi 'l. 


4.3 Influence of the cosmological model 

The probability for a source with a known apparent flux and 
redshift to be gravitationally lensed with the formation of 
multiple images depends on the cosmological model. This 
may be seen for instance through the dependence of the in¬ 
finitesimal light-distance element cdt/dz on the cosmological 
mass density parameter in Eq. as well as through the 
dependence of the lensing cross section on (see Eqs. 
and 0 via the definition of the angular diameter distances. 
If we consider a flat expanding FLRW universe model, the 
only dependence of the lensing probability on the universe 
model is made through the cosmological mass density pa¬ 
rameter Q.m- Indeed, although a flat expanding FLRW is to- 



Figure 6. Average lensing probability as a function of the cos¬ 
mological matter density parameter Q,m- 


tally characterised by Q,m and Hq, an increase in the value of 
Hq will only act as a scaling factor (decreasing the lensing 
volume while increasing the density of deflectors). So far, 
we have considered a FLRW fiat universe with Q.m = 0.3 
and Ho = 70 km s“^ Mpc“^. In Fig. [6l we display the be¬ 
haviour of the average lensing probability (r) for the XXL 
sources as a function of fim, for a flat universe with Ho = 70 
km s“^ Mpc“^. We have computed both (tsie) and (tsis)- 
Both models lead to the same behaviour of the mean lensing 
probability as a function of Qrn and, here as well, we observe 
that the SIS mass distribution leads to a slight overestimate 
of the average lensing probability when compared to that 
corresponding to the SIE model, whatever the value of flm- 
There is a very strong dependence of the expected 
fraction of lensed sources in the survey on Qrn- For this 
reason, the statistics of gravitational lensing in a well- 
defined sample of sources has been widely u sed to probe 


1984 iFukueita et al.lll990j. 

lTurneilll99d. 

Keeton 19981. Ghae et al.l 

I 2 OO 2 . Keeton 


The XXL sample on its own will not allow better con¬ 
straints on the value of flm than re cent lens surveys such as 
the SDSS-LQS. IOguri et al.l ll2012lf for instance constrained 
^Lambda tO Q.Lambda = 0.79^ ^ ± 0.06(syst.) OU the 

basis of 19 lenses from the SDSS-QLS statistical sample. 
The XXL lensed source sample will thus be combined with 
other recent surveys (including the SDSS-QLS) to better 
constrain the cosmological parameters. 


5 CONSIDERING THE R-BAND CUT-OFF 

In the previous chapter, when calculating the amplification 
bias in Eq.0 we have only considered the distribution of the 
sources as a function of their flux fx in the X-ray band and 
we have considered the lensing optical depth r ( 2 s, fx) to be 
only a function of the source redshift and its flux in the X-ray 
band. Ideally, the source r-magnitude in the optical domain 
should also be considered for the calculation of r (2s, fx, r) 
and dobs (za, fx,r) and the calculation of the mean lensing 
optical depth (r) in Eq. should include the integretation 
over r. 

The joint probability density dobs {za, fx, r) may be de¬ 
composed as 

doba [Za,fx,r) = doba {Za , fx) d {r\Za , fx) (15) 

where d {r\za, fx) is the normalised distribution as a func¬ 
tion of r for sources with a redshift and an X-ray flux 
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Figure 7. Left : distribution of the COSMOS sources in the (/x?’") plane. We have represented the X-ray fluxes in terms of the X-ray 
magnitudes rrix = —2.51og/x, and have represented the X-ray and r-band cut-off of the XXL survey. We have also represented the 
contour plot of the COSMOS source density in the {fx^p) plane. Right : amplification bias as a function of the total amplification of 
the lensing event when considering only the DNCF as a function the X-ray flux fx (continuous dark grey line) and when averaging over 
the r-band magnitude distribution of the sources. 


in the ranges [zs, + dzs] and [fx, fx + dfx], respectively, 
and Eq. HD remains correct if we consider the optical depth 
{zs, fx) averaged over the r magnitudes 

'^{r) (zajx) = J T{zs,fx,r)d{r\zs,fx)dr. (16) 

Equivalently, (zs, fx) can be calculated through Eq.|3]if 
we consider an average lensing cross section 

S{r> JJ B^^)dy, (17) 

Sy 


where we have defined the amplification bias averaged 
over the r-magnitudes 


B{r) — 


/ 


N {zs, fx/Ax,r -\-2.5log{Ar)) 
N{zs,fx,r) 


dir\zs,fx)dr, 


(18) 

where N {zs, fx, r) represents the density of sources with a 
redshit Zs , having an observed flux fx in the X-ray and a r- 
band magnitude r, Ax and Ar are the amplifications due to 
the lensing event in the X-ray and r-band, respectively. If the 
source is point-like, we may then assume that Ax = .4^. If 
not, the ratio Ax jAr will depend on the size of regions in the 
AGN emitting the X-ray and the optical fluxes, respectively, 
as well as the positions of these regions with respect to the 
caustic curves. The expectation of the ratio Ax jA-r will thus 
necessitate heavy simulations to be performed. As a first 
approximation we consider point-like sources. 

In this work, the amplification bias Bx-ray has been 
calculated so far uniquely on the basis of the X-ray flux 
through the relation 


Bx 


— ray 


Njfx/Ax) 

N{fx) 


(19) 


where we have neglected the redshift dependence of the 
sources as a function of their redshift because of a too small 
source sample to correctly characterise this possible redshift 


dependence. In the following we will compare the amplifica¬ 
tion bias obtained through Eqs. [18] and 1191 

To calculate rigorously the amplification bias through 
Eq. [18] we thus need to determine the distribution of the 
source population in the (za, fx,r) space. Our observational 
data from the COSMOS sample is not large enough to esti¬ 
mate this distribution over the whole area probed by the sur¬ 
vey. Consequently, we assume the distributions N (za, fx, r) 
and d{r\za, fx) to be independent of the redshift. 

Let us now estimate the distribution N {fx,r) on the 
basis of the distribution of the COSMOS sample in the 
{fx,r) plane, represented on the left panel of Fig. [ 7 ] We 
have represented all the COSMOS sources for which the 
soft X -ray flux and r-ban d magnitude is given in the sam¬ 
ple of iBrusa et al.l 1 I 2 OIOI '). As 98% of the detected X-ray 
sources have an optical counterpart in the r-band for X-ray 
sources brighter than the XXL 0.5 detection probability cut, 
we can consider this sample to be a complete sample of X- 
ray sources brighter than the X-ray cut, for which we know 
the optical counterpart. 

To estimate N {fx,r) we have proceeded as follows. We 
bin in terms of the = —2.5 log fx magnitude, with a bin 
width of drria: = 0.5 and bin centres ranging from 33.5 to 37. 
In each rrix bin, we construct the histogram of the sources as 
a function of the r-band magnitude, using bins with a width 
of dr = 0.5 and bin centres ranging from 17 to 28 with step 
of 0.5. Each histogram as a function of r is fitted by means 
of a Gaussian profile. The number of sources in each bin is 
then divided by drrixdr. 

We thus obtain the evolution of the Gaussian fit param¬ 
eters (i.e. the amplitude Ac, the average position re and the 
standard deviation ac) as a function of the bin and we 
fit the dependence of these parameters as a function of rrix 
by a linear law in the rrix range 33.5 to 37. We model the 
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density of the sources in the {f^, r) plane by 
N Ux,r) = Ao (m.) exp 

The COSMOS data do not have enongh bright sonrces to 
constrain N {fx, r) for brighter than 33.5. We have rep¬ 
resented the iso-density contonrs of the calculated function 
N {fx,r) on the left panel of Fig. [71 showing a very good 
agreement between the observed distribution of the COS¬ 
MOS sources and the modelled density function. 

We now use the modelled function N {fx, r) to calculate 
the amplification bias in Eq. 1181 for sources with different 
X-ray to optical flux ratios. Assuming the amplification in¬ 
duced by the lensing event in the X-ray and the optical are 
identical, when amplified, a source is displaced in the {rrix, r) 
parallel to the dashed line shown on the left panel of Fig. (T) 
Consequently, the amplification bias is calculated thanks to 
the evolution of N (fx, r) along this trajectory. 

On the right panel, we represent the behaviour of the 
amplification bias as a function of the total amplification of 
the lensing event, for 4 sources with apparent X-ray fluxes 
equal to the XXL X-ray limiting flux fxxL-cut, and 5, 10 
and 15 times the value of the fxxL-cut- The dashed light 
grey curve represents the amplification bias B x-ray ob¬ 
tained when only considering the DNCF of the XXL-sources 
as a function the X-ray flux, calculated through Eq. 1191 The 
continuous black curve corresponds to the average bias 
calculated thanks to Eq. 1181 We see a perfect agreement be¬ 
tween these curves and the amplification bias calculated only 
considering the DNCF as a function of the X-ray flux fx- 
Consequently, the amplification bias calculated through Eq. 
[ 6 ] perfectly corresponds to the amplification bias averaged 
over the r-band magnitudes. 

Thus, assuming a point-like source, we may calculate 
the amplification bias in the combined X-ray/optical data 
by considering uniquely the X-ray distribution of a complete 
and deeper sample, which validates the method introduced 
in Section (2] Nevertheless, in the analysis of the final XXL 
sample, we will have to also consider the redshift dependence 
of the amplification bias, which will be made possible thanks 
to the much larger size of the sample. 


6 CONCLUSIONS 

We have calculated the expected statistical properties of the 
multiply imaged sources to be detected among the optical 
counterparts of the XXL point-like sources, modelling the 
deflectors successively with SIS and SIE mass profiles. We 
find that 

• one expects the formation of 20 (21 using the SIS 
model) multiply-imaged AGNs out of which 16 (17 for the 
SIS case) should be detected among the optical counterparts 
with an angular resolution of 0.45”, and we only expect the 
detection of one gravitational lens system composed of more 
than 2 lensed images; 

• the expectations are consistent when modelling deflec¬ 
tors with SIE and SIS mass distributions, although the SIS 
model leads to a slight overestimate of the mean lensing 
probability. This overestimate is a function of the amplifica¬ 
tion bias and is thus different for each survey; 


• the late-type galaxy population should not contribute 
to the lensed sources to be detected; 

• although the detection is done simultaneously in the X- 
ray and in the optical domain, the amplification bias may 
be estimated from the X-ray flux distribution, as long as we 
consider a complete X-ray sample from a deeper survey and 
for point-like sources. 

In this work, we have considered isothermal profiles to 
model the deflectors. This has allowed us to get first good 
estimates of the expected number of lensed AGN in the XXL 
survey. However, more detailed calcu lation s ought to be car¬ 
ried o ut. Indeed. [Auger et"^ (|2010l 'l and iKoopmans et al.l 
ll2009ll . through the analysis of massive early-type deflec¬ 
tors from the SLAGS survey, have found a slight deviation 
from the isothermal profile, with a steeper slope parame- 
ter (( 7 ) = 2.078 and 2.085, respectively, where the mass 
distribution evolves as r ’'). If we consider two mass dis¬ 
tributions with a same total mass, a steeper profile would 
lead to a higher fraction of the lens in the center, which will 
increase the Einstein angular radius, therefore leading to an 
increase in the lensing cross section of the deflector. This in¬ 
crease in the lensing cross section with the steepness of the 
radia l mass profile has been pu t in evidence in a series of pa - 
pers (iMandelbaum et al.ll2009l and Ivan de Ven et ahllioo^ i. 
where the authors studied the impact of galaxy shape and 
density profile on the selection biases in surveys for the de¬ 
tection of strong lenses. Because isothermal profiles are sin¬ 
gular the authors analysed more realistic profiles in order to 
define the total m ass. _ 

Furthermore, ISonnenfeld et all l|2013l ') have shown from 
the study of the SL2S Galaxy-scale Lens Sample (with de¬ 
flectors in the range 0.2 < 2 : < 0 . 8 ) that the mass density 
profile of early-type galaxies depends on the redshift, lower 
galaxies showing a steeper average profile. As deflectors with 
a steeper profile tend to have a higher lensing cross section, 
this would favour the deflector redshift distribution to be 
shifted towards lower redshift. 

These effects, along with the redshift dependence of the 
amplification bias will have to be considered in the analysis 
of the final XXL sample. 
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